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1.  Introduction 


All  materials,  due  to  their  atomic  nature,  have  some  magnetic  response.  However, 
for  most  materials,  the  response  is  so  small  that  it  is  considered  negligible.  Magnetic 
susceptibility,  %,  is  a  material  property  that  indicates  how  strongly  a  material  will 
magnetize  in  response  to  an  applied  magnetic  field.1  Consequently,  it  is  also  a 
measure  of  the  field  strength  required  for  producing  an  appreciable  response.  For 
materials  with  exceptionally  small  magnetic  susceptibilities,  a  significantly  strong 
magnetic  field  is  necessary  to  produce  a  response  (2-10  T).2-6  This  is  the  case  for 
many  ceramic  materials  (e.g.,  alumina),  which  are  often  paramagnetic  or 
diamagnetic.  2-6 

Many  ceramics  exhibit  magnetic  susceptibility  anisotropy  (A/)  with  variations  in 
magnetic  susceptibility  due  to  different  crystallographic  directions.  Because  of  this 
anisotropy,  crystals  will  align  their  axes  with  the  highest  magnetic  susceptibility  in 
the  direction  of  the  field  that  provides  a  mechanism  for  preferred  orientation.  3  9  The 
magnetic  response  of  ceramics  is  governed  by  2  main  equations  derived  from  first 
principles  according  to  Sugiyama  et  al.10  Equation  1  describes  the  anisotropic 
magnetic  energy  (AE),  or  torque,  acting  on  the  particle  due  to  its  magnetic 
anisotropy  as  a  function  of  the  system  properties  and  applied  magnetic  field.10 
Equation  2  describes  the  temporal  requirement  ( t )  for  rotation  of  the  crystal, 
assuming  no  steric  hindrance.10  It  is  derived  by  summing  the  effects  of  inertial, 
viscous,  Lorentzian,  and  magnetization  forces.  Equation  2  assumes  that  the 
characteristic  inertial  response  time  is  negligible  compared  with  the  viscous  and 
Lorentzian  times,  which  is  valid  for  most  small  particles  with  moderate  densities 
and  viscosities.10  Full  derivations  are  available  in  Sugiyama  et  al.,10  and  an 
abbreviated  derivation  and  adaptation  of  relevant  equations  is  available  in  the 
Appendix  for  Eqs.  1  and  2. 

AE  =  V  (*c~*a-b)g2  >kBT  (1) 

2Ho  B  W 


and 


t  =  “ho 


30ri+r2aB2  . 

5 (Xc~Xa,b)B2  U 


(  tand 


\tan6 


;) 


(2) 


where  V  is  the  primary  crystal  volume,  X“,b  and  x<-  are  the  unitless  magnetic 
susceptibilities  of  the  crystallographic  planes  that  are  perpendicular  and  parallel  to 
the  magnetic  field  following  alignment  respectively,  pio  is  the  permeability  of  free 
space,  B  is  the  externally  imposed  magnetic  flux  density,  kn  is  Boltzmann’s 
constant,  T  is  the  absolute  temperature,  t  is  the  time  needed  for  a  crystal  to  rotate  a 
certain  amount  assuming  no  steric  hindrances,  r)  is  the  viscosity  of  the  melt  or  fluid 
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surrounding  the  crystal,  r  is  the  crystallite  radius,  o  is  the  electrical  conductivity  of 
the  particle,  and  0o  and  0  are  the  initial  and  final  angles  of  the  crystal  easy-axis 
(parallel  axis,  x<)  with  respect  to  the  magnetic  field,  respectively.10  Historically, 
there  is  much  confusion  in  the  literature  concerning  these  equations  due  to  the  use 
of  several  different  sets  of  units  for  magnetic  properties.  The  equations  and 
derivations  referenced  in  this  report  were  confirmed  for  accuracy  by  performing 
dimensional  analysis. 

Using  these  equations,  the  magnetic  torque  applied  to  the  particle  and  the  time 
necessary  for  rotation  can  be  estimated  and  used  to  predict  the  response  of  materials 
to  strong  magnetic  fields.  The  sensitivity  of  such  responses  to  small  changes  in 
input  variables  such  as  field  strength  and  viscosity  can  also  be  assessed.  This  report 
presents  predictions  about  the  magnetic  field  response  and  sensitivity  of  process 
variables  for  alumina  ceramics  that  are  currently  being  investigated.  However,  the 
code  and  analysis  developed  for  this  application  can  easily  be  adapted  to  other 
material  systems  by  changing  the  appropriate  material  properties. 

2.  Methods 


2.1  Defining  Parameters 

For  each  process  variable,  representative  values  were  gathered  experimentally, 
from  the  literature,  and  from  product  data  sheets.  Ay  for  pure  alumina  was 
calculated  from  the  molar  susceptibility,  Aym,  which  was  found  to  be  7.1  x  10'8 
emu/mol.5  Since  emu/mol  is  a  nonstandard  centimeter-gram-second 
electromagnetic  unit  (emu),  it  was  first  converted  to  SI  (international  standard) 
units  of  m3/mol  by  multiplying  by  the  conversion  factor  471  x  10"6.11  Finally,  molar 
susceptibility,  Aym,  was  converted  to  unitless  volume  susceptibility,  Ay,  using 
Eq.  3. 12 

AXm  =  AXgFW  =  JFW  ^  Ax  =  =  VmAxm  ,  (3) 

where  FW  is  formula  weight,  p  is  the  theoretical  density  of  the  material,  Vm  is  the 
molar  volume  of  the  material,  and  Ayg  is  the  mass  susceptibility.12  Values  for 
alumina  were  set  at  0.10196  kg/mol  [2(26.98)  +  3(16.00)  g/mol]  and  3990  kg/m3, 
respectively.13 

Particle  radius,  r,  was  set  at  150  nm  to  represent  the  nanosized  powder  currently 
being  investigated.  Particle  volume,  V,  was  calculated  from  the  particle  radius 
assuming  spherical  geometry.  Temperature,  T,  was  set  at  room  temperature 
(298  K)  and  electrical  conductivity,  a,  was  set  to  0  (ohm-m)"1  since  the  conductivity 
of  alumina  is  negligible.14  Applied  magnetic  field  strength,  B ,  and  medium 
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viscosity,  rj,  were  analyzed  at  2  different  values  based  on  practical  testing 
conditions.  Field  strength  was  analyzed  at  both  1.8  and  9  T  (magnetic  fields 
available  for  testing),  and  viscosity  was  analyzed  at  both  1.215  and  0.75  Pa-s16  to 
mimic  casting  with  epoxies  and  gel-casting  slurries  as  the  liquid  mediums, 
respectively.  0n  and  0  were  set  to  89°  and  1°,  respectively.  Considering  that  the 
range  of  possible  angles  before  symmetry  causes  repetition  (in  terms  of  energy)  is 
0°-90°  (see  Fig.  1),  it  was  determined  that  the  maximum  alignment  time  will  occur 
when  all  particles  are  at  90°  to  the  field  direction.  However,  since  the  tangent 
function  is  undefined  at  90°,  the  initial  angle  was  set  to  89°  so  that  useful  values 
could  be  calculated.10  For  similar  reasons,  the  final  angle  was  set  to  1°  instead  of 
0°.10 


Fig.  1  Schematic  displaying  the  convergence  of  angles  with  respect  to  the  magnetic  field  (0) 
due  to  symmetry 

2.2  Magnetic  Response  Predictions 

The  infonnation  and  equations  presented  were  coded  into  MATLAB  functions  for 
easy  evaluation.  Viscosity,  field  strength,  particle  radius,  and  magnetic 
susceptibility  anisotropy  (associated  with  changes  in  doping  levels)  were  varied 
while  holding  all  other  variables  constant  within  ranges  relevant  to  the  conditions 
presented  for  alumina  with  respect  to  alignment  time. 

2.3  Sensitivity  Analysis 

Using  the  same  MATLAB  functions  from  Section  2.2,  sensitivity  analysis  was 
perfonned.  Sensitivity  analysis  is  commonly  used  when  evaluating  mathematical 
models  of  complex  biological  enzymatic  pathways,  but  the  principle  is  still  useful 
for  this  application.17,18  Sensitivity,  in  modeling  terms,  is  defined  as  the  relative 
change  in  output  compared  with  the  relative  change  in  input.  Higher  values  of 
sensitivity  are  associated  with  larger  impact  for  a  specified  variable.  While  this 
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indicates  a  lack  of  robustness  in  biological  models,  for  mathematical  models 
derived  from  first  principles  it  identifies  the  most  efficient  parameters  to  target 
experimentally.18  In  this  analysis,  each  variable  was  systematically  perturbed  by 
decreasing  it  by  10%  while  holding  all  other  variables  constant.  Since  sensitivity 
analysis  can  only  be  conducted  locally,  relative  to  a  baseline  set  of  parameters, 
sensitivity  was  assessed  at  both  1.8  and  9  T  and  0.75  and  1.2  (N-s)/m2  to  determine 
if  the  local  sensitivity  differed  experimentally  from  purely  theoretical  inputs  at 
values  of  interest.  The  times  required  for  alignment  or  magnetic  torque  energy  were 
recorded  both  before  and  after  the  perturbation,  and  the  sensitivity  function,  SOF, 
was  calculated  according  to  Eq.  4. 

CQp  (OV-tput Perturbed- OntPUt Baseline) 

( OutputBaseune )  x(perturbation  amount )  ’  v  ' 

where  output  indicates  either  AE  (magnetic  torque  energy)  or  t  (alignment  time), 
and  perturbation  amount  is  the  amount  that  the  variable  is  reduced,  which  in  this 
case  is  0.1  (10%)  for  all  analyses.17  A  value  of  10%  perturbation  was  used  to 
maintain  locally  relevant  conditions  within  the  model,  as  required  when  performing 
sensitivity  analysis.18 

3.  Results  and  Discussion 


3.1  Magnetic  Response  Predictions 

The  effects  of  viscosity  on  alignment  time  with  varying  applied  magnetic  field 
(Fig.  2)  and  the  changes  in  magnetic  anisotropy  with  varying  magnetic  field 
(Fig.  3)  have  been  studied.  As  magnetic  field  increases,  alignment  occurs  more 
rapidly.  Increasing  viscosity  also  increases  the  time  necessary  for  alignment  but  is 
not  as  effective  as  the  influence  of  a  magnetic  field.  Increasing  magnetic  anisotropy 
also  decreases  the  alignment  time  significantly,  as  expected.  However,  since  it  is 
unknown  how  much  the  magnetic  anisotropy  of  alumina  is  affected  by  different 
dopants  and  dopant  concentrations,  it  is  also  unknown  how  relevant  this  will  factor 
will  be  experimentally. 
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Effect  of  Viscosity  on  Alignment  Time 


Fig.  2  Effect  of  viscosity  (ij)  on  alignment  time  with  varying  applied  magnetic  field  ( B ).  For 
relevant  system  conditions  of  a  1.8-T  applied  magnetic  field  and  a  viscosity  of  0.75-1.2 
(N-s)/m2,  time  necessary  for  alignment  is  predicted  to  be  in  the  range  of  35-55  s. 


Fig.  3  Effect  of  magnetic  anisotropy  (A/)  on  alignment  time  with  varying  applied  magnetic 
field  ( B ).  For  relevant  system  conditions  of  a  1.8-T  applied  magnetic  field  and  a  Ax  of  8.92  x 
10'12,  time  necessary  for  alignment  is  predicted  to  be  in  the  range  of  50-60  s. 
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Figure  4  shows  the  effect  of  particle  radius  on  alignment  time  with  varying 
electrical  conductivity.  As  conductivity  increases,  alignment  time  increases  with  an 
increase  in  the  Lorentzian  force  component  that  opposes  rotation.  The  Lorentz  force 
arises  when  current  induced  by  the  magnetic  field  interacts  with  the  field  itself, 
inducing  electromagnetic  force  on  the  particles.  This  force  works  in  the  opposite 
direction  of  the  magnetic  torque  forces  and  can  be  thought  of  as  a  “magnetic  inertia” 
that  must  be  overcome  for  alignment  to  occur.  The  induced  current  increases  with 
increasing  electrical  conductivity  such  that  higher-conductivity  materials 
experience  higher  Lorentz  forces.  However,  it  should  be  noted  that  when  the 
conductivity  is  zero,  as  in  the  case  of  alumina,  particle  radius  has  no  effect  on 
alignment  time. 


Effect  of  Particle  Radius  on  Alignment  Time 


r  (nm) 

Fig.  4  Effect  of  particle  radius  (r)  on  alignment  time  with  varying  particle  conductivity  (<r). 
For  relevant  system  conditions  where  particle  conductivity  is  0  or  nearly  0,  particle  radius  is 
predicted  to  have  no  effect. 

While  it  is  unknown  how  much  dopant  identity  and  concentration  will  affect  overall 
conductivity,  the  values  presented  in  Fig.  4  are  well  outside  the  range  of  any 
reasonable  change,  particularly  for  most  ceramic  materials  of  any  composition.19 
The  reason  such  a  large  variety  of  conductivity  values  was  used  was  to  produce 
changes  in  alignment  time  significant  enough  to  show  the  effect  of  particle  radius 
and  conductivity.  Looking  at  the  scale  on  the  y-axis  in  Fig.  4,  alignment  time  only 
varies  by  12  x  10'11  s  with  a  100-fold  change  in  conductivity.  This  indicates  that  in 
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the  case  of  almost  all  ceramic  materials,  the  Lorentzian  force  component  of  the 
magnetic  alignment  equations  (and  associated  variables)  can  be  considered 
negligible  as  has  been  proposed  previously.20  This  is  particularly  evident  in  the 
context  of  parameters  such  as  magnetic  field  strength  and  anisotropy,  which  have  a 
much  stronger  influence  over  alignment  time.  Additionally,  simulations  were 
performed  for  both  silver  and  graphene,  which  are  known  as  some  of  the  most 
conductive  materials  available  (a  =  3.5  x  107  [ohm-m]"1  and  2.05  x  107  [ohm-m]'1, 
respectively21),  and  changes  in  alignment  time  due  to  particle  radius  were  still  only 
on  the  order  of  2  s.  For  this  reason,  it  was  concluded  that  for  powder  samples  of 
almost  any  material,  variations  due  to  changes  in  conductivity  and  particle  radius 
can  be  neglected,  especially  when  the  set  times  are  on  the  order  of  tens  of  minutes. 

For  the  levels  of  applied  magnetic  field,  viscosity,  magnetic  anisotropy, 
conductivity,  and  particle  radius  relevant  to  alignment  of  alumina  in  a  1.8-T 
magnetic  field  with  low-viscosity  slurries,  the  simulations  predict  necessary 
alignment  times  of  approximately  30-50  s,  which  is  well  within  the  range  of 
experimental  conditions.  However,  the  simulations  do  not  take  into  account 
possible  steric  interactions  between  particles  that  have  been  shown  to  prevent 
complete  alignment  for  irregularly  shaped  particles.22  This  is  an  issue  particularly 
when  using  particles  of  plate  or  needle-like  morphologies  that  can  pin  each  other 
during  the  rotation  process,  creating  particle  bridges  that  prevent  further  rotation. 
This  inhibitory  effect  will  only  increase  with  the  solids  loading  of  the  slurries,  as 
there  will  be  even  less  room  for  particles  to  freely  rotate. 

In  addition,  the  simulations  do  not  address  the  possibility  that  agglomerates  have 
many  crystals  with  opposing  domains  that  may  conflict  with  each  other  during 
alignment.  Since  the  current  powders  are  not  milled  or  de-agglomerated  in  any 
significant  way  before  processing,  it  is  very  likely  that  they  are  highly 
agglomerated,  resulting  in  less-effective  alignment.  These  simulations  also  do  not 
take  into  account  changes  in  magnetic  torque  based  on  variations  in  particle  shape, 
which  is  known  to  affect  magnetic  response.23  25  Given  that  the  predicted  times  are 
on  the  order  of  tens  of  seconds,  it  is  still  anticipated  that  some  degree  of  alignment 
will  be  achievable,  as  the  setting  times  for  epoxy  and  gel-casting  slurries  range  from 
30  min  to  several  hours.  This  is  also  supported  by  preliminary  results.26 

3.2  Sensitivity  Analysis 

The  calculated  sensitivity  functions  (SOFs)  for  parameters  in  the  magnetic  torque 
energy  and  alignment  time  equations  (Eqs.  1  and  2)  are  shown  in  Tables  1  and  2 
for  the  1.8-  and  9-T  magnetic  fields  and  the  0.75-  and  1.2-(N-s)/m2  baseline 
viscosities,  respectively.  In  both  models,  the  local  sensitivity  is  exactly  the  same 
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for  all  cases  of  experimental  interest.  This  was  expected  because  the  model  was 
derived  from  first  principles  instead  of  empirically,  which  can  introduce 
unexpected  and  unrealistic  outputs  in  extrapolated  regions.  Instantaneous 
sensitivity  (not  nonnalized  to  the  baseline)  was  calculated  analytically  at  values  of 
1.8  T  and  0.75  (N-s)/m2  using  partial  derivatives.  The  resulting  SOF  values  (data 
not  shown)  showed  similar  trends  but  much  higher  magnitudes  compared  with  the 
data  in  Table  2,  which  supports  the  validity  of  this  method.  However,  since  the  local 
nonnalized  sensitivity  values  relative  to  variations  that  could  occur  naturally  in  a 
practical  setting  are  much  more  valuable  than  instantaneous  theoretical 
sensitivities,  the  method  previously  presented  is  considered  to  be  more  relevant. 

Table  1  Sensitivity  factors  for  process  variables  in  the  magnetic  torque  energy  equation 
(Eq.  1).  Higher  values  indicate  higher  model  sensitivity,  which  highlights  variables  that  are 
more  effective  for  increasing  material  response  to  a  magnetic  field.  Ax  denotes  the  magnetic 
susceptibility  anisotropy,  r  denotes  particle  radius,  and  B  denotes  applied  magnetic  field 
strength. 


Variable 

Sensitivity  Function  (SOF) 
1.8  T  9  T 

1 

1 

r 

2.71 

2.71 

B 

1.9 

1.9 

Table  2  Sensitivity  factors  for  process  variables  in  the  alignment  time  equation  (Eq.  2). 
Higher  values  indicate  higher  model  sensitivity,  which  highlight  variables  that  are  more 
effective  for  increasing  the  material  response  to  a  magnetic  field.  B  denotes  applied  magnetic 
field  strength,  //  denotes  slurry  viscosity,  r  denotes  particle  radius,  and  Ax  denotes  the  magnetic 
susceptibility  anisotropy. 


Variable 

Sensitivity  Function  (SOF) 

1.8  T  9  T 

0.75  (N-s)/m2 

1.2  (N-s)/m2 

0.75  (N-s)/m2 

1.2  (N-s)/m2 

B 

2.35 

2.35 

2.35 

2.35 

> 1 

1.00 

1.00 

1.00 

1.00 

r 

0.00 

0.00 

0.00 

0.00 

1.11 

1.11 

1.11 

1.11 

As  shown  in  Table  2,  when  considering  the  parameters  that  affect  alignment  time, 
the  strength  of  the  applied  magnetic  field  ( B )  has  the  most  influence.  Magnetic 
anisotropy  (A/)  and  viscosity  (//)  also  affect  the  alignment  time  but  not  as  drastically 
as  field  strength.  This  was  also  reflected  in  Figs.  1  and  2,  as  described.  The  particle 
radius  has  a  sensitivity  factor  of  zero  and  therefore  has  no  effect  on  alignment  time 
when  particles  are  nonconductive  as  discussed  previously. 

While  of  little  experimental  relevance,  the  sensitivity  of  magnetic  torque  energy  to 
various  process  variables  is  displayed  in  Table  1.  In  contrast  with  alignment  time, 
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particle  radius  has  the  highest  impact  on  magnetic  torque,  followed  by  field  strength 
and  magnetic  anisotropy.  While  seemingly  counterintuitive,  the  mathematical 
reason  for  this  phenomena  is  easily  explained.  When  volume  is  expressed  as  a 
function  of  r3  and  combined  with  the  viscous  and  Lorentz  force  terms,  it  is  simply 
canceled  out,  as  shown  in  the  Appendix  (Eqs.  10  and  1 1).  Conceptually,  the  reason 
is  not  as  obvious.  Volume  does  not  translate  directly  to  shorter  alignment  times 
because  its  advantage  is  offset  by  a  viscous  force  term  in  which  the  resistance  to 
torque  is  proportional  to  the  increase  in  volume.  Therefore,  this  term  does  not 
appear  in  the  final  alignment  time  equation  (Eq.  2). 

4.  Conclusions 


Simulations  predict  that  alignment  will  be  possible  under  the  imposed  experimental 
conditions.  Furthermore,  when  dopants  are  added  to  enhance  the  magnetic 
anisotropy  (Ax),  the  alignment  response  will  increase.  Previous  literature  cites  field 
strengths  in  excess  of  10  T  necessary  for  alignment,5,6  but  according  to  these 
simulations,  it  appears  that  lower  strength  fields  are  sufficient  for  alignment. 
Preliminary  experimental  evidence26  and  some  literature27  appear  to  support  this,23 
but  more  work  is  required  to  confirm  that  lower  fields  are  indeed  sufficient  to 
observe  statistically  significant  alignment  for  alumina-based  ceramics.  It  is  possible 
that  the  equations  derived  from  first  principles  simply  do  not  take  into  account  other 
practical  considerations,  such  as  the  aforementioned  interparticle  interactions,22  and 
conclusions  regarding  alignment  time  are  not  applicable. 

The  sensitivity  analysis  indicates  that  the  strongest  parameter  affecting  magnetic 
alignment  is  the  strength  of  the  applied  magnetic  field  ( B ).  This  means  that  in 
experimental  applications  the  most  effective  way  to  increase  alignment  is  to 
increase  the  magnetic  field.  However,  this  also  means  that  every  precaution  must 
be  taken  to  ensure  that  field  is  stable  and  consistent  before  commencing 
experiments.  At  the  very  least,  the  instability  of  the  field  must  be  monitored  to 
quantify  uncertainty,  as  this  is  likely  the  most  prevalent  source  of  variability. 
Viscosity  and  magnetic  susceptibility  anisotropy  (A/)  were  also  shown  to  have  an 
impact  on  overall  alignment,  but  their  effects  were  less  pronounced  than  changes 
in  field  strength.  While  it  was  not  a  source  of  variability  within  experiments, 
analysis  of  the  sensitivity  indicated  that  changes  in  anisotropy  affect  alignment 
response,  thus  variations  in  alignment  at  lower  field  strengths  with  different  dopants 
are  also  expected. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


J{a,b,  (Ka.b) 


Xc  (Kc) 


Ax  (Ak) 


B 

T 

>1 

t 

k-B 

jUo 

r 

a 

V 


0 


0o 


AE 

AXm 

AXg 

Vm 

P 

FW 

SOF 


volume  magnetic  susceptibility  of  crystallographic  planes 
perpendicular  to  the  applied  magnetic  field 

volume  magnetic  susceptibility  of  crystallographic  plane 
parallel  with  the  applied  magnetic  field 

volume  magnetic  susceptibility  anisotropy  associated  with  a 
material,  calculated  as  Xc  -  Xa.b, 

applied  magnetic  field  [T] 

temperature  [K] 

solution  viscosity  [(N-s)/m2] 

time  needed  for  crystal  rotation  [s] 

Boltzmann’s  constant,  1.381  x  1 0"23  [J/k] 

penneability  of  vacuum,  47t  x  10"7  [(T-m)/A] 

crystal  radius  (assumes  spherical  particles)  [m] 

electrical  conductivity  of  the  particles  [(D-m)"1] 

primary  crystal  volume  [m3] 

final  angle  of  magnetic  easy-axis  (c-axis)  with  respect  to  magnetic 
field  axis  [°] 

initial  angle  of  magnetic  easy-axis  (c-axis)  with  respect  to 
magnetic  field  axis  [°] 

energy/torque  associated  with  magnetic  anisotropy  [J]  or  [N-m] 

molar  magnetic  susceptibility  anisotropy  associated  with  a 
material  [m3/mol] 

mass  magnetic  susceptibility  anisotropy  associated  with  a 
material  [m3/kg] 

molar  density  [mol/m3] 

theoretical  density  associated  with  a  material  [kg/m3] 

formula  weight  associated  with  a  material  [kg/mol] 

sensitivity  function:  how  much  that  parameter  affects  the 
outcome  of  a  model 
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Appendix.  Deviation  of  Magnetic  Alignment 
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As  adapted  from  Sugiyama  et  al.,1  4  kinds  of  forces  are  active  when  aligning 
particles  in  a  magnetic  field:  1)  inertial  force,  2)  viscous  force,  3)  Lorenz  force,  and 
4)  magnetization  force.  Equations  for  each  of  these  active  forces  are  presented  in 
Eqs.  A-l  through  A-8. 


Inertia  Force-. 

2  r  d2e 

sPr  7^- 

(A-l) 

Viscous  Force-. 

8nrjr3  ^ . 

(A-2) 

Lorenz  Force-. 

4  r  dd 

— nr 5  —  . 

15  dt 

(A-3) 

Magnetization  Force-. 

—  VAxB2sin(20). 

2  Ho 

(A-4) 

Summing  all  torque  forces  acting  on  particles  under  a  magnetic  field  results  in 
Eq.  A-5,  which  is  consequently  solved  to  produce  the  equation  for  alignment  time 
as  presented  in  Eq.  2  in  the  main  report.  Eq.  A-5  is  simply  the  sum  of  main  report 
Eqs.  5-8. 

-pr5  ^  +  8nrir3  —  +  —  nr3  —  +  —r3AxB2  sin(20)  =  0.  (A-5) 

5  t/  t  t/ 1,  t/ 1  2  |f  Q 

As  stated,  since  the  particles  are  small  and  the  viscosity  and  particle  material 

density  in  question  (alumina)  moderate,  the  inertia  force  term  can  be  neglected 

since  under  those  circumstances  the  characteristic  response  time  will  be  much 
smaller  than  the  contribution  of  the  other  forces.  Neglecting  the  inertial  term  yields 
Eq.  A-6: 

8nrir3  —  +  —nr3  —  +  —  r3AxB2  sin(20)  =  0  ,  (A-6) 

dt  15  dt  2  \Xq 

which  upon  simplification  results  in  Eq.  A-7,  which  is  the  instance  where  particle 
volume,  expressed  as  r3  in  the  magnetization  term,  is  mathematically  removed  from 
the  alignment  time  equation. 

8n^n+TS,Tr2n  +  ir/xB2sm(2e)  =  0.  (A-7) 

Solving  the  differential  equation  yields  Eq.  A-8,  which  when  the  appropriate 
substitutions  are  made,  exactly  corresponds  to  Eq.  19  in  Sugiyama  et  al.,1  where 
the  authors  have  chosen  to  lump  all  material  parameters  into  a  tenn,  r. 

tan(0)  =  tan(f>„)  exp  (-  ^ 

Eq.  A-8,  upon  rearranging  to  solve  for  t,  yields  Eq.  2  in  the  main  report. 


'Sugiyama  T,  Tahashi  M,  Sassa  K,  Asai  S.  The  control  of  crystal  orientation  in  non-magnetic  metals  by 
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